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on the ionization of molecules by a strong, linearly polarized laser pulse. Electron angular distri- 
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the molecular tunneling theory appropriately extended to non-linear polyatomic molecules. Our 
investigations show that the orientational dependence of ionization yields is primarily determined 
by the nodal surface structure of the molecular orbitals. 
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I. INTRODUCTION 

Ionization is an ubiquitous process in the interaction 
between atoms and molecules and intense laser fields. 
The ejected electron propagates in the combined laser 
and Coulomb field, and the laser field may steer back the 
electron and force it to recoUidc with the parent ion. In 
this way, the release of an electron into the laser-dressed 
continuum can initiate strong-field processes, such as 
above threshold ionization, multiple ionization, and high 
harmonic generation (see, e.g., Ref. Q and references 
therein). Clearly, the laser field only drives the elec- 
tron in the direction of the laser polarization vector 
and, consequently, rescattering dynamics is typically only 
present for linearly polarized laser fields. The possibility 
of rescattering depends not only on the parameters of 
the laser. Also the scattering angle with which the elec- 
tron leaves the parent ion plays a role. For most atomic 
ground states, strong-field ionization leads to an outgoing 
electronic current in the direction of the linear polariza- 
tion of the laser, and hence a recoUision is possible after 
the field has changed its direction. For molecules, how- 
ever, it happens, that the molecular symmetry imposes 
one or more nodal planes in the ground state wave func- 
tion. In such planes the ejection of photoelectrons may 
be strongly suppressed Q , causing recoUisionally induced 
double ionization not to peak for linearly but for ellipti- 
cally polarized laser light It is therefore clear that a 
good understanding of the dependence of the initial ion- 
ization process upon laser parameters and upon param- 
eters characteristic for the system, such as symmetry, is 
crucial for the further development and possible practical 
utilization of strong-field phenomena. 

A large number of studies has established that the 
laser parameters determining strong-field ionization are 
the wavelength, the intensity, and the polarization P, 0, 
0, 0, 0- We note in passing that in the case of few- 



cycle laser pulses, the ionization rate and the photoelec- 
tron angular distribution also depends critically on the 
phase of the carrier frequency with respect to the en- 
velope j^. For atoms the most important system pa- 
rameter is the binding energy of the initial state. In 
fact, the simple ADK tunneling model 9j, employing 
only the laser intensity and the atomic binding energy 
as input parameters, provides for many atomic species a 
very good description of measured total ionization rates. 
Also, the theor y re ferred to as the strong-field approx- 
imation (SFA) lli m m m m has been successful 
in modeling atomic ionization including the angular dis- 
tribution of the ejected photoelectrons. In this theory, 
the ionization takes place from a field-free ground state 
to a laser-dressed continuum final state via the interac- 
tion with the lowest-order laser-atom operator. In the 
final state, the Coulomb interaction between the electron 
and the parent ion is neglected and the electron is de- 
scribed by the state of a free electron in the laser field - 
the so-called Volkov state (see Ref. for a very recent 
discussion). 

For molecules a few recent studies have shown that 
there is an additional system parameter that also strongly 
influences strong-field ionization namely the s ymm etry of 
the initial electronic state 0,111 E 111 Ellll IH M 

. Modifications of both the ADK tunneling theory |23 
and the SFA 0| have been introduced to account for the 
molecular symmetry. Hereby the two models have pro- 
vided results in reasonably agreement with measured ion- 
ization rates of molecules. In the molecular ADK (MO- 
ADK) theory '2^ one expands the asymptotic molecular 
ground state wave function in a one-center basis and ap- 
plies the atomic tunneling theory to each partial wave. 
The atomic SFA is straightforwardly generalized to the 
molecular SFA (MO-SFA) by using molecular wave func- 
tions. A very appealing feature of MO-SFA is that the 
expansion of the molecular wave function in linear combi- 
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nation of atomic orbitals (LCAO), allows an interpreta- 
tion of the spectrum of photoelectrons in terms of inter- 
ference from the multi-center molecule; an interpretation 
which in the context of singic-photon ionization was of- 
fered almost 40 years ago . 

In our previous work Q we used the molecular strong- 
field approximation (MO-SFA) in the velocity gauge for- 
mulation to illustrate the strong dependence of the pho- 
toelectron angular distribution and the total ionization 
rates of ethylene and benzene on (i) the orientation be- 
tween the molecule and the laser polarization and on (ii) 
the symmetry of the electronic wave function. The pur- 
pose of the present paper is to give a more thorough dis- 
cussion of the influence of molecular symmetry on strong- 
field ionization. In particular, we study four different 
molecules representing four different classes of symmetry 
(point groups). We show how a reduction of the molec- 
ular symmetry lifts the degeneracy of molecular orbitals 
and leads to different ionization channels with distinct 
features. 

More recently, we found that the MO-SFA in the ve- 
locity gauge formulation failed to explain orientational 
effects in ionization of N2. However, using the length- 
gauge to describe the molecule-field interaction, results in 
good agreement with experiments 0| were obtained . 
Consequently, in the present work we found it necessary 
to employ both gauges and check if the results obtained 
are consistent. In addition, we calculate ionization rates 
by the MO-ADK tunneling model and compare with the 
results from the MO-SFA. 

The paper is organized as follows. In Sec.|nl we present 
a justification for using a Hartree-Fock calculation to rep- 
resent the electronic ground states of the molecules under 
study. In Sec. UTTI we recall the theory of the MO-SFA 
and the MO-ADK tunneling models. In particular, we 
present the generalization of the MO-ADK theory |2j| to 
the case of non-linear polyatomic molecules. In Sec. IIVI 
we provide details of the calculations and discuss the 
choice of basis functions, the choice of nuclear coordi- 
nates, and the method for obtaining the ionization rates. 
In Sec. E] we present the numerical results comprising a 
comparison between the length and the velocity gauge. 
We argue that the length gauge gives the better results 
and use it together with the MO-ADK model to calculate 
ionization rates for each of the four molecules. In Sec. lVll 
we conclude and provide an outlook for the future. 



while the HOMO inferred from experiments is a ag or- 
bital jl^]. When including electron correlations in the 
calculations, the correct ordering and ionization energies 
are obtained [2^. Another example where the effects of 
electron-electron correlation are important is F2 which 
is not even bound in a HF treatment. It is therefore 
not surprising that one-electron models have failed in 
the attempt of explaining strong-field ionization of F2 

diiiiia. 

In this work, we consider ionization of ethylene and 
substituted benzenes. For these systems electron cor- 
relation may be safely neglected in the description of 
the initial state, since photoionization spectra in the va- 
lence region of these molecules consist of clearly sepa- 
rated bands which are unambiguously assigned and are 
in accordance with Hartree-Fock predictions .22J . In 
these systems, each electron occupies an orbital obtained 
from a HF calculation and it is an accurate approxima- 
tion to describe ionization by the removal of an electron 
from an occupied orbital while leaving the ion in an un- 
relaxed state. With these approximations, the many- 
electron problem translates into the much simpler single- 
electron problem. Only ionization from the highest or 
next highest occupied molecular orbital (referred to as 
HOMO-1) is considered. For ethylene and benzene the 
separations between the first and second ionization po- 
tential are 2.3 eV and 2.2 eV, respectively. In Ref. 
it was shown that such large energy differences led to 
a much lower ionization probability from the HOMO-1 
than from the HOMO. The separation between the first 
and second ionization potential is, however, much lower 
for the substituted benzenes, approximately 0.4 eV. Ac- 
cordingly, we will consider ionization from the HOMO 
only for ethylene and the degenerate HOMO of benzene 
while for the substituted benzenes we will additionally 
include ionization from the HOMO-1. Ionization from 
all other orbitals is neglected. 

III. THEORY 

In this section, we describe the theories that we apply 
for the calculation of the molecular strong-field ioniza- 
tion. Throughout this work, we fix the nuclei at their 
equilibrium positions as determined with the HF theory. 



II. ELECTRONIC STRUCTURE: 
HARTREE-FOCK OR BEYOND 

In some molecules the inclusion of electron correlations 
is necessary in order to predict even qualitatively correct 
properties. For example, in N2 the ordering of the or- 
bitals obtained from a simple Hartree-Fock (HF) calcu- 
lation is not consistent with the observed single-photon 
photoelectron spectrum: The calculated highest occu- 
pied molecular orbital (HOMO) of N2 is a 7r„ orbital 



A. Molecular strong-field approximation 

The molecular strong-field approximation (MO-SFA) 
0, 0, llOl is a generalization of the atomic strong-field 
approximation 0, ITTl IT^ . For a linearly polarized elec- 
tric field F{t) = Fq sin (uit) with associated vector po- 
tential A{t) — Aq cos (ust), Aq = Fq/oj, the laser electron 
interaction is represented by [atomic units {h — nif, = 
e = ao = 1) are used throughout] 

y^i^«)(i) = r.F(t), (1) 
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in the length gauge (LG), and by 



(2) 



in the velocity gauge (VG). 

In either gauge, we utilize the 27r/ijj = T periodicity of 
the field to express the angular differential, dW/dq, and 
total, W, ionization rates as sums of n-photon absorp- 
tions 13, 



dW 
dq 



2tT \Aqnfqn 



W ^ 2^Y1 / \Aqn\''qndq, 



(3) 
(4) 



l/T {■^v{r,t)\V^'\t)\^{r,t))dt (c ^ {LG,VG}) are 
evaluated at the momentum g„ = ^J2{nuj — E}, — Up), 
with Eh the binding energy of the initially bound 
electron and Up the ponderomotive potential. Energy 
conservation also determines the minimum number of 
photons needed to reach the continuum, i.e., the lower 
limit uq of the sum. In the expression for Aq„, \E'y is a 
(27r)~'^/^ normalized Volkov wave function and ^'(r,i) 
the wave function for the electron in the combined field of 
the molecule and the laser. Introducing the strong-field 
approximation, i.e., 5'(r,t) ~ \l/o(r',t) = e'^''*$o(^), 
with $0 the solution of the time-independent field-free 
HF equation, we find the length and velocity gauge 
amplitudes pollpl 



where the transition amplitudes Aq^i 



/((LG) ^ ± 

qn fj-i 



dt -Eh - 



Qlit) 



^o{Qn{t))expi ncut 



qn ■ OiQ sm(ujt) + — sm(2ujt) 
2lo 



(5) 



/l(VG) ^ 
qn 



dt exp i nujt 



qn ■ a.Q siniut) + — sm(2ujt) 
2lu 



(6) 



Here Qn{t) = qn + A{t) denotes the time-dependent 
momentum, cxq — Aq/uj denotes the quiver radius, and 
$g(q) = (27r)^'^/^ J dre~^'^"^^o{r) the momentum wave 
function of the initially bound electron. The velocity 
gauge transition amplitude of Eq. © may be simplified 
by replacing the time integral with a generalized Bessel 
function Jn{u,v) |l2j | and using energy conservation 



transform of the wave function is then 

MQu) = E E <^^J4>^J(.qn)e~^''-''\ (9) 

i=l j=l 

with representing 4>ij in momentum space. 



Ag""^ ^{Up-niu)^o{qn)J-r 



qn ■ "0, 



Up 
2uj 



(7) 



Applying Eq. iQ to the differential rate of Eq. © is 
equivalent to the expression used in Ref. 17J. We note 
that Eq. ((Sj simplifies to Eq. ^ when neglecting A{t) 
compared to q„, i.e., when Qn(t) « q„. 

We construct the molecular one-electron wave func- 
tions as linear combinations of basis functions ipij based 
on each atom (LC AO-MO), located at the positions Ri. 
The wave function can now be written as 



N J 



Mr, {R^}) =E Yl'^Ajir - R,), (8) 

i=i j=i 

where i runs over the N atoms in the molecule and j runs 



over the basis functions on each atom 



The Fourier 



B. Molecular tunneling theory 

The molecular ADK (MO-ADK) tunneling theory 
[19I I23II is a generalization of the atomic ADK theory 
[9|. In previous works [H 13, El the MO-ADK theory 
was applied to diatomic molecules. Here the theory is 
generalized to cover also nonlinear polyatomic molecules. 

The tunneling theory relies on the assumption that at 
any given instants of time the molecule will respond to 
the external laser field as if it were a static electric field. 
The rate of ionization in the oscillating field may thus be 
determined by the time-averaged static rates. Whether 
this approach is reasonable or not depends on the value 
of the Keldysh parameter 7 = \J2Ebjr |3 ^'^^^ 7^1 
in the tunneling regime. 

The formulation of a tunneling theory in the case of 
molecules is complicated compared to the atomic case by 
two related features. Firstly, the presence of multiple nu- 
clei within the molecule breaks the spherical symmetry 
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of the field-free system and necessitates a description of 
the molecular wave function in terms of a superposition 
of partial waves. Secondly, the Euler angles describing 
the orientation between the laboratory fixed frame (with 
a Z-axis determined by the linear polarization vector of 
the external field) and the molecular body fixed frame 
have to be specified. Note that we use {X, Y, Z) to la- 
bel the laboratory fixed frame of reference and {x, y, z) 
to label the body-fixed frame of reference. The tunnel- 
ing rate of molecules can be determined once the field- 
free asymptotic wave function is known. In a body-fixed 
frame, labeled by superscript B, this function must follow 
the asymptotic Coulomb form 



where 



Qmllm(r), (10) 



I,' 



where Zion is the charge of the residual ion and k is re- 
lated to the binding energy, k — \/2Eb. For the tun- 
neling process in the DC case, we assume the electric 
field to point in the positive Z-direction, corresponding 
to a situation where tunneling would occur in the nega- 
tive Z-direction. Consequently, we need to express the 
asymptotic form of the molecular wave function in that 
direction. If the body-fixed frame is rotated by the Euler 
angles ((/), 9, x) with respect to the laboratory fixed sys- 
tem, the asymptotic form in the laboratory fixed system 
of the field-free molecular ground state wave function 
is expressed through the rotation operator D{(f),9,x) as 
^o('') — -0(0, 0, x)$^(r), where we have used the su- 
perscript L to denote the laboratory fixed system. Equa- 
tion (|I0|) then leads to the expression 



l.m 



^(0 



(11) 

where ^m'm('/', ^, x) is a Wigner rotation function '53, 
Is^ . In Eq. ((TT|l the sum over m' and the correspond- 
ing Wigner functions describes the rotation between the 
coordinate systems and the sum over Z,m is a signature 
of the breaking of the spherical symmetry by the molec- 
ular system. For linear molecules the projection m of 
the electronic angular momentum onto the body-fixed 
axis is a good quantum number and, hence, for such sys- 
tems there would be no summation over m in Eqs. H1U|) 
and Hll|l (23] . For later convenience, we note that </> and 
X represent rotations around the space-fixed Z-axis and 
the body fixed z-axis, respectively, while 9 represents the 
angle between the Z and z axes. 

From the asymptotic form of Eq. , the total ioniza- 
tion rate in a static (DC) field in the positive Z direction 
is calculated as in the atomic case 0) HE > and the 
result is 



(12) 



2k' 



3 s 2Z/K-\m'\-l 
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In a slowly varying field, the ionization rate is found 
by averaging the DC rate over an optical cycle 



W = 



1 

2^ 



2tt 



Wstat (Fq COs(wt)) d{LUt). 



(14) 



The DC rate is given by Eq. (|12|l when the field is oriented 
in the positive Z direction, |7r — cj^j > 7r/2, corresponding 
to tunneling in the negative Z-direction. When the field 
points in the negative Z direction, 7r/2 < ojt < 37r/2, cor- 
responding to the possibility for tunneling in the positive 
Z-direction, the geometry is equivalent to a field pointing 
in the positive Z direction but with an inverted molecule. 
We account for this by applying the parity operator 



on the wave function 



(r) 



l.m 



G™(-l)'r/,„(r), (15) 



and we see by comparing with Eq. IjlOl) that this simply 
corresponds to the substitution Cim (— l)'C;m in Eq. 
()13|l . Note that parity eigenstates, corresponding to in- 
version symmetric molecules, contain only even or odd I 
states. When the field direction is changed Eq. (|13|l will 
either remain invariant or change sign and the DC rate is 
thereby left invariant. Contrary, the DC rate will not be 
invariant to field inversion for states which are not parity 
eigenstates. This is also to be expected since the wave 
functions in the tunneling regions in the positive and neg- 
ative Z direction will be different. Under the assumption 
/ Fq ^ 1 the integral in Eq. (|14|l may be approximated 

by 



W = 



3FoW+M) + W-,,,iFo) 



(16) 



where Wg^^^^^Fo) are the DC rates for the positive and 
negative field directions with respect to the Z direction. 



IV. CALCULATIONS 

A. Basis set 

The Hartree-Fock wave functions for ethylene, ben- 
zene, fluorobenzene and o-chlorofluorobenzene are calcu- 
lated using GAMESS [s^ with a standard valence triple 
zeta basis set |3 . This basis set contains s and p orbitals 
for carbon, fluorine and chlorine while the hydrogen ba- 
sis consists only of s orbitals. Since the tunneling rates 
and the length gauge MO-SFA rates both rely on the 



asymptotic form of the wave function a precise descrip- 
tion is needed in this region. This is obtained by adding 
an extra diffuse s and p basis function. The orbitals of 
the molecules considered are odd under reflection in their 
molecular planes. Thus the LCAO expansion, Eq. ((SJ, 
can only include p orbitals orthogonal to the plane. For 
example, if we define the yz plane to coincide with the 
molecular plane, the relevant p orbitals are the Px or- 
bitals. 



B. Expansion and asymptotic coefficients 

The calculation of the initial molecular ground state 
wave function is performed at the equilibrium geome- 
try predicted by Hartree-Fock theory. Bond lengths ob- 
tained in this way are typically within few pm from ex- 
perimental data [221 ■ The nuclear geometries are shown 
in Fig. [2 

Note that we define the body-fixed z axis to be per- 
pendicular to the molecular plane for benzene while it 
lies in the molecular plane for the remaining molecules. 
Since we want to investigate the influence of molecu- 
lar alignment on the ionization rates, we have chosen 
our axes to coincide with the principal axes for each 
molecule. Aligned molecules are generated as superpo- 
sitions of rotational eigenstates which are most conve- 
niently described in this coordinate system ^3 ■ Our def- 
initions of coordinates are in accordance with the usual 
conventions except for chlorofluorobenzene. 

The asymptotic form of the orbitals, Eq. (|10|) . is found 
by expanding in spherical harmonics 

<fo(r-)=^^^i™(r)yi™(f), (17) 

l,m 

leading to 

Fi„,{r) = I dnYi:Ar)Mr). (18) 

These functions are then matched to the form Fim[r) ^ 
Cimr^'°°^"'~^e~'^'^ in the asymptotic region whereby the 
coefficients Cim are determined. The values of Cim for 
the molecules investigated in this work are given in Ta- 
ble HI All the orbitals are odd eigenstates of the re- 
flection operation and they are additionally chosen purely 
real. The projection on the spherical harmonics, Eq. (|18|l . 
will then be either purely real (m odd) or purely imagi- 
nary (m even) . Another consequence of the yz reflection 
antisymmetry is Fi^-m{f) = —Fi^mif)- In ethylene, ben- 
zene and fluorobenzene our orbitals are additionally be 
eigenstates of the Uxz reflection operation. If the wave 
function is an even (odd) eigenstate of the Gxz reflection 
operation then only the real, m odd, (imaginary, m even) 
radial functions will occur. 

Since the basis functions we use are Gaussian and the 
correct asymptotic form to which we fit, is exponential, 
we cannot expect a precise fit all the way to infinity. 
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Figure 1: Nuclear geometries and isocontours for the active 
orbitals of (a) ethylene, (b) benzene, (c) and (d) HOMO 
and HOMO-1 for fluorobenzene respectively and (e) and 
(f) HOMO and HOMO-1 for 1,2-chlorofluorobenzene respec- 
tively. The origins of the coordinate systems are in the center- 
of-mass. For later convenience we see that all the orbitals 
contain nodal planes in the molecular plane. The HOMOs of 
fluorobenzene and 1,2-chlorofluorobenzene, (d) and (f), con- 
tain nodal surfaces which are approximately parallel to the xy 
plane. The HOMO-1 of fluorobenzene contains a nodal plane 
in the xz plane while the HOMO-1 of 1,2-chlorofluorobenzene 
contains a nodal suface which is almost in the xz plane. The 
relation between the body-fixed axes shown in the figures and 
the laboratory-fixed axes is discussed in Sec. 1111 51 



Accordingly, we have adjusted the range of asymptotic 
fitting to start well outside the nuclear positions and is 
limited such that the radial functions Fim (r) of Eg. H18|) 
still follow the exponential behavior. In Ref. [23, we 
obtained accurate results by this procedure. 



C. Simplifications in the length gauge 

In the length gauge MO-SFA the transition amplitude, 
will only depend on the asymptotic behavior of the coor- 
dinate space wave function [23, |31| ■ 

When we have determined the angular coefficients Cim 
we may simply use the Fourier transform of Eq. Hll|) in 
the length gauge transition amplitude of Eq. jSJ. The 
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Table I: Asymptotic coefficients Cj, 



/ 


m 


C2H4 


CsHo 


C6H5F 


1,2 C6H4CIF 










HOMO 


HOMO-1 


HOMO 


HOMO-1 


1 


±1 


±1.10 




±0.49 




±1.27 


±0.08 


2 


±2 








Tl.70i 


T0-39j 


Tl.39i 




±1 




±1.36 


Tl.70 




=F0.95 


±0.19 


3 


±3 






±0.13 




=F0.28 


=F0.54 




±2 








±1.22i 


±0.39j 


±1.5H 




±1 


±0.22 




±0.89 




±1.83 


=F0.71 


4 


±4 








±0.35j 


±0.102 


±0.4H 




±3 






±0.24 




±0.13 


±0.16 




±2 








T0-68i 


T0.63i 


T0.75i 




±1 




=F0.32 


±0.31 




=F0.53 


±0.26 


5 


±5 










±0.03 


±0.11 




±4 








T0-25i 


TO.IOj 


TO.SOi 




±3 






=F0.11 




=F0.21 


=fo.oi 




±2 








±0.33i 


±0.062 


±0.30i 




±1 






±0.23 




±0.82 


=F0.25 


6 


±4 








±Q.16j 


±0.10i 


±0.18i 




±2 








T0.12i 


T0.22i 


TO.lOi 




±1 










=F0.17 


±0.07 



Fourier transform of Eq. is 



Ix)Yi„Aq), (19) 



where the radial momentum-space functions fi{q) are 

r(/±f + 2) 



(20) 



X 2i^l 



Z±#±2;±#±3, 3 



with 2^i(a, fe; c; z) Gauss' hypergeometric function |4lj| . 
Equation (jSJ can now be rewritten as 



l.7n,m' 



(21) 



X - / dti-Et 



)Il{Qn{t))yinAQnit)) 



orientations (0, 9, x) by multiplication with the appro- 
priate Wigner functions. We have checked this method 
against calculations with the fully numeric Fourier trans- 
forms applied in Eq. Q and we find agreement within 
< 20% in the final rates. 

In closing this section, we note that with the above sim- 
plifications of the length gauge MO-SFA, the only non- 
trivial dependency on electronic structure of the MO- 
SFA LG and the MO-ADK theory is through the dm 
expansion coefficients. 



D. Ion signal 

The ionization rates in Eqs. Q and (|14|) depend on 
the instantaneous amplitude of the field and the molec- 
ular orientation described by the Euler angles ((/>, 6*, x)- 
We only consider linearly polarized light and hence the 
results are independent of </>, the angle of rotation around 
the polarization vector. For a Gaussian laser beam with a 
Gaussian temporal profile with FWHM r, the amplitude 
of the vector potential is 

A(7^,Z,0 . ^^e-V.(.)%,p 

(22) 

where {TZ, Z) are the cylindrical coordinates, Iq the peak 

intensity, wq is the spot size and w{Z) = itig \/l~+^Z^7^5 
where Zfj = ttwq/X is the Rayleigh length and A the 
wavelength. 

The rate equations for the ionization probability of a 
molecule oriented according to the Euler angles (6*, x) and 
located at (TZ, Z) in the laser focus are 

^ = -P^(A(7^,Z,^),0,x)po (23) 

^ = w{A{n,z,t),e,x)po. (24) 

where po and pi denote the probabilities of having a neu- 
tral or an ion, and ionization to higher charge states is ne- 
glected. By the end of the pulse, the solution to Eq. 
is given by 

Pl(7^, Z, 9, x) = 1-exp j W (A(7^, Z, t'), 9, x) dt' 

(25) 

The orientational dependent number of ionized molecules 
A^i is found by integrating Eq. (|25|) over the beam profile 



exp ^ (^nu;t + q„ • ao sin(c.i) ± g sin(2..t)) , N,ie,x) ^ 2np J UdU J dZp, {n,Z,9,x). 



(26) 



with Qn{t) = Qn + A{t). This formula is particularly use- 
ful when we want to consider many different molecular 
orientations, since the integrals in Eq. 1)21(1 are indepen- 
dent of molecular orientation. Firstly, we need to calcu- 
late all the integrals once for each I and m' at some lab 
fixed momentum q„. Then we can easily get the transi- 
tion amplitude at the same lab fixed momentum for all 



where p is the constant density of the target gas. In ex- 
periments it is difficult to measure the absolute yield due 
to unknown detection efficiency. The measured ion sig- 
nal is, however, proportional to the number of ionized 
molecules. Measured ratios of yields for different molec- 
ular orientations are thus independent of detection effi- 
ciency. We note that the actual evaluation of the integral 
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Eq. (|25|l can be simplified by integrating over isointensity 
shells 42]. This method describes experiments in which 
all ionized molecules are measured. 



V. RESULTS AND DISCUSSION 

To justify the use of the approximate theories discussed 
in Sec. IIIII we first calculate the intensity dependent 
ion signal for ethylene and benzene and compare with 
experimental data. Previous studies using the velocity 
gauge MO-SFA have shown that the absence of the fi- 
nal state Coulomb interaction in this model will under- 
estimate the ionization rates. To take this interaction 
into account a correction factor C^^^j — (k^/Fo)^'^'""/'* 
was proposed In our velocity gauge calculations, we 
adopt this correction factor. A recent study on N2 j2flj re- 
vealed that the effect of the Coulomb interaction is much 
smaller with the length gauge MO-SFA. The explanation 
is that the combined laser-ion interaction with the elec- 
tron is already present through the factor $o('Zn + ^(0) 
in Eq. © and also in the length gauge, the ionization 
occurs at large distances from the core. In this spatial re- 
gion, the Coulomb interaction is already suppressed com- 
pared to the electron-laser interaction. Below we demon- 
strate that the length gauge MO-SFA without any extra 
Coulomb correction factor reproduces experimental ion 
signals of C2II4 and CeHg. 

For direct comparison of results obtained by MO-SFA 
and MO-ADK theories with the experimental data it is 
necessary to express the total rate of ionization as the 
sum of rates leading to unresolved final states. This 
means that the rates in Eqs. Q and must be multi- 
plied by the number of electrons which occupy the active 
orbital, i.e., 4 for the degenerate HOMO of benzene and 
2 for the non-degenerate closed shell orbitals. 

Figures [3 and 13 show theoretical results and experi- 
mental data of the ion signal versus peak intensity for 
experiments on ethylene and benzene, respectively. The 
calculations are the MO-ADK, the length gauge MO-SFA 
(MO-SFA LG) and the velocity gauge MO-S FA (MO-SFA 
VG). Following the discussion in Sec. lIVDl we normalize 
the results of each calculation to match the experimental 
data at high intensities. The figures show that the MO- 
SFA LG predicts results in excellent agreement with both 
the ethylene and the benzene experimental data. The 
Coulomb corrected MO-SFA VG results fit the data rea- 
sonably well, while we find a poor agreement between the 
MO-ADK model and experiments. We explain this lat- 
ter discrepancy by the fact that the intensities considered 
here do not strictly correspond to the tunneling regime. 
Typical values of the Keldysh parameter are 3.0 (2.8) for 
ethylene (benzene) at /q = lO^^ Wcm'^ and 0.94 (0.88) 
at /o = lO^Wcm"^ We note that despite both MO- 
SFA models are in agreement with these total ion-yield 
experiments, a recent study on orientational-dependent 
ionization from N2 showed that the length gauge formu- 
lation is the better choice [20l |. 




Peak intensity (Wcm^^) 



Figure 2: Ion signal for ethylene at 800 nm versus the peak 
intensity of the laser. The laser wavelength is 800 nm, the 
beam waist is 49 \rai and the pulse duration (FWHM) is 200 fs. 
Experimental data are from Ref. [i^l . The solid line indicates 
the the length gauge MO-SFA (MO-SFA LG) predictions, the 
long-dashed line shows the velocity gauge MO-SFA (MO-SFA 
VG) predictions and the short-dashed line shows predictions 
of MO-ADK theory. 




Peak intensity (Wcm"^) 



Figure 3: Similar to Fig. |5| but for benzene. Experimental 
data are from Ref. |4^. 



A. Ethylene 

The main conclusion from our earlier velocity gauge 
work on ethylene was that the angular distributions 
of photo electrons are very sensitive to the relative orien- 
tation between the laser polarization and the molecule. 
This could be interpreted as an effect of the nodal plane 
structure of the HOMO. Electron emission was found to 
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Figure 4: (a)-(d) Polar plots of the angular differential ioniza- 
tion rate of ethylene by an 800 nm, 5 x 10^'' Wcm~^ field us- 
ing the length gauge MO-SFA. Rates are in units of 10^" s~^. 
The polarization (double headed arrow) is horizontal and the 
molecular plane is perpendicular to the paper with the C-C 
axis (dotted line) in the plane of the paper {x = 0°) and (a) 
6» = 0°, (b) 30°, (c) 60°, (d) 90°. 



be impossible in the directions of the molecular plane 
which led to a suppression of the total ionization rate 
when the molecular plane coincided with the polariza- 
tion vector. 

We first present polar plots of angular differential ion- 
ization rates using the length gauge MO-SFA (Fig.^ and 
the Coulomb corrected velocity gauge MO-SFA (Fig. [SJ 
for different molecular orientations, specified by the an- 
gle 9 between the polarization and the body-fixed z axis 
defined in Fig. ^ (a) . Our primary purpose of these cal- 
culations is to check if the results in the two gauges are 
compatible. From Figs.21and|Slwe see that both models 
predict an increase in the ionization rate as the polar an- 
gle 9 increases. When the molecular plane coincides with 
the polarization, the rates are strongly suppressed and in 
both models no electron emission is observed along the 
polarization direction. In the velocity gauge this phe- 
nomenon was explained by the presence of a nodal plane 
in the molecular plane There is a one-to-one corre- 
spondence between nodal planes in coordinate space and 
momentum space. Thus, the velocity gauge transition 
amplitude, Eq. lO, vanishes at momenta which lie in the 
nodal plane. This is indeed what we observe in Fig. 
In the length gauge transition amplitude, Eq. ||SJ), the 
momentum wave function is not evaluated at the actual 
outgoing momentum but at a time-dependent displaced 
momentum q„ -|- A{t). This has the consequence that 
even though q„ lies in the nodal plane the transition am- 
plitude can be nonzero. Only when q„ and A{t) both lie 
in the nodal plane {9 = 0°) the transition amplitude is 
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Figure 5: Similar to Fig. 0] but calculated in velocity gauge. 



zero. This argument is general for molecules which ion- 
ize from orbitals containing nodal planes. In both gauges 
the ionization rate is thus expected to be suppressed for 
molecular orientations with the polarization coinciding 
with a nodal plane. This leads to the qualitative agree- 
ment in the two gauges. For molecules without planar 
nodal surfaces no such suppressions can a priori be ex- 
pected for any geometry. Based on previous work on 
strong-field ionization of N2 0| and detachment of neg- 
ative ions 31, 45] we judge that the length gauge is the 
proper choice of gauge to use in the strong-field approx- 
imation. Accordingly, in the calculations below we will 
only use the length gauge MO-SFA. 

Currently, angle resolved photoelectron measurements 
on aligned molecules have not been performed. How- 
ever, the first measurements of total ion yields on aligned 
molecules have been reported 16]. In order to make pre- 
diction on these types of experiments we show in Fig. El 
the total ion yield as a function of molecular orientation 
using (a) the MO-SFA LG model and (b) the MO-ADK 
model. The figure shows the yields which would be mea- 
sured if all molecules were oriented according to the Eulcr 
angles specified. Explicitly, we show the signals as a func- 
tion of the polar angle 9 for two specific x angles with 
all signals being normalized to the 9 = 0° geometry. We 
recall that 9 = 0° for the molecular plane parallel with 
the linear polarization, and 9 = 90° for the perpendicu- 
lar geometry. In experiments, utilizing linearly polarized 
light pulses for the alignment, the molecular axis with 
largest pola rizability is aligned along the polarization di- 
rection [ig (the C-C axis for ethylene). In this case the 
molecule is free to rotate around its body-fixed z axis 
and measured signals are thus a result of an average over 
the X angle. Correspondingly, we have calculated the x 
averaged signal too. 
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Figure 6: The orientational dependent ion signal for ionization 
from ethylene using (a) the molecular strong-field approxima- 
tion in the length gauge formulation and (b) the molecular 
tunneling theory. The signal is given as a function of the po- 
lar angle with x = 0° (long dashed), x = 90° (short dashed) 
and X averaged (solid). The parameters of the laser are: 
A = 800 nm, Jo = 5 x 10^^ Wcm"^, r = 20 fs and wq = 49 fim. 



Figure El(a) shows, in accordance with Fig. 01 a large 
increase in the ion yield with increasing polar angle for 
X = 0° as well as for the x averaged result. For x = 90° 
the polarization vector is always directed in the molec- 
ular plane and by the argument given above the yield 
is generally much lower. In Fig. (b) we see much the 
same trend for the MO-ADK model with the x = 90° 
geometries being suppressed compared with their x = 0° 
counterparts. Contrary to the results from the MO-SFA 
LG model, the yield is not maximized when polarization 
is perpendicular to the plane for x = 0° and x aver- 
aged. The theories disagree because the tunneling rate is 
derived under the assumption that the electron emission 
only occurs close to the space-fixed Z axis determined by 
the linear polarization direction of the field. Hence, the 
MO-ADK rates reflect only the situation in this spatial 
region. For atoms, this is a reasonable assumption since 
the initial spherical symmetric potential will be affected 
most strongly in this direction and the potential barrier 
which arises from the combined field and atomic potential 
will then be minimized. For molecules, the initial poten- 
tial is more complicated and the potential barrier can be 
minimized in other directions. Important contributions 
to the ion yield are thereby not taken into account and 
this is the reason why the results of the theories do not 
agree. 

Since the tunneling ionization occurs close to the Z 
axis by assumption, the interpretation of the MO-ADK 
results can be made by analyzing the electron distribu- 
tion in this direction. Previously, it was noted that the 
tunneling rate is largest when the initial electronic cloud 
is aligned with the field direction [S^- In Fig. [7| (a) we 
show a contour plot of the HOMO in the body fixed xz 
plane. When the field vector lies in this plane the geom- 
etry corresponds to x = 0° and with vertical and hori- 
zontal directions corresponding to 6* = 0° and 9 = 90°, 
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Figure 7: (a) Contour plot of the HOMO of ethylene in the xz 
plane, (b) The norm square of the HOMO along the directions 
e = 50° and e = 90°. The unit of length is the Bohr radius 
ao = 5.29177 x 10"" m. 



respectively. To make it clear that 9 = 90° does not 
correspond to the largest electron density we have fur- 
thermore made a plot of |<I>o(r)p along the directions 
9 = 50° and 9 = 90° (Fig.[71(b)). We clearly see that the 
electron density is largest for 9 = 50° at large distances. 
The minimum around 9 — 90° is also present in Fig. |^ 
(b) for the x averaged yields. 



B. Benzene 

Benzene is a symmetric top molecule and hence it can 
only be aligned in the molecular plane |46j while it will 
be free to rotate around the body-fixed z axis (see Fig. Q] 
(b)). Thus only x averaged signals as shown in Fig.|Hlcan 
be measured. We choose to present our calculated sig- 
nals as a function of the angle, O, between the molecular 
plane and the polarization. Note that this angle is related 
to the usual Euler angle through the relation O = 90° — 9 
where 9 is the polar angle between the Z-axis and the 
moleular plane in the frame defined in Fig.^ This choice 
of data representation is convenient since a comparison 
with benzene derivatives like fluorobenzene and chloroflu- 
orobenzene can be made more directly (Sees. IV (]l and 
IVDj. In Fig. m we see that the MO-SFA LG model pre- 
dicts the maximum yields at = 42° and O = 138° 
while the MO-ADK model predicts maximum yields at 
O = 33° and O = 147°. In both models the yield is min- 
imized at = 0°, 90° and 180°. These findings are in 
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Figure 8: The orientational dependent ion signal for ionization 
of benzene as a function of = 90° — 6. The parameters of 
the laser are given in the caption of Fig. [5] 



accordance with our expectations since for 8 = 0°, 
90° and 180° the polarization vector lies in a nodal plane 
and the asymptotic electron density is maximized in be- 
tween. In a quantitative comparison between the models 
we find a somewhat larger modulation in the orienta- 
tional dependent yield with the MO-ADK model than in 
the MO-SFA LG model. 



Fluorobenzene 
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When one hydrogen atom in benzene is replaced by a 
fluorine atom, the six-fold rotational symmetry is broken 
and consequently the degeneracy of the HOMO is lifted. 
The doubly degenerate HOMO of benzene sphts up in 
two components with the ionization potentials 9.35 eV 
and 9.75 eV. This energy difference is small compared 
with the energy separations to the other occupied orbitals 
and we will thus consider ionization from both orbitals. 

Using the MO-SFA LG model, we present in Fig.|ni(a)- 
(c) the ion yields which originate from (a) the HOMO, (b) 
the second highest occupied orbital (H0M0~1) and (c) 
the integral yield from both channels. The figures|ni(d)- 
(f) contain equivalent results from the MO-ADK model. 
We find that the signals from both orbitals contribute 
significantly to the total signal and depending on the ge- 
ometry, the contribution from either one will dominate. 
The relative importance between the ionization channels 
at a given geometry can be understood as a consequence 
of the symmetries of the orbitals. 

Due to the C2v symmetry of fluorobenzene the orbitals 
are eigenstates of the xz reflection operator with eigen- 
values -t-1 (even) or —1 (odd). The reflection symmetry 
of an orbital is indicated by the subscript 1 for the even 
states and 2 for the odd states. The HOMO of fluo- 
robenzene is a hi orbital. Fig. ^(c), while the HOMO-1 
is an a2 orbital. Fig. ^(d). Both orbitals are odd with 
respect to reflections in the molecular {yz) plane. Simi- 
lar to ethylene, ionization of fluorobenzene is suppressed 
when X = 90° compared with the x average since the 



Figure 9: The orientational dependent ion signal for ioniza- 
tion from the highest occupied orbitals of fluorobenzene using 
MO-SFA LG (a)-(c) and MO-ADK (d)-(f). See Fig.Elfor de- 
tails on the laser pulse and legends. 



molecular plane is a nodal plane for both active orbitals. 
The bi orbital contains a non-planar nodal surface which 
is approximately parallel to the xy plane. When 6 = 90° 
the polarization is accordingly nearly parallel to this sur- 
face and we find the yield from the HOMO to be mini- 
mized around = 90°. The 02 orbital is odd under the 
xz reflection and thus it contains a nodal plane in the 
xz plane. When the polarization is parallel to this plane 
(x = 0°) we also find suppression of the ionization sig- 
nal. In all cases we find qualitative agreement between 
the MO-SFA LG and MO-ADK results. In order to dis- 
cuss the effect of the substitution of hydrogen with flu- 
orine we note that the angles O for benzene and 9 with 
X = 0° for fluorobenzene both describe the angle between 
the respective molecular planes and the polarization. It 
is thus relevant to compare the curve in Fig. |S1 with the 
dashed curve in Fig. |51 In both theories we find that the 
total yield for fluorobenzene is slightly more suppressed 
around 9 = 90° than for benzene when it is compared 
with the maximum yields around 45 degrees (MO-SFA 
LG) and 30 degrees (MO-ADK). 
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Figure 10: Similar to Fig.^lbut for chlorofluorobenzene. 



D. Chlorofluorobenzene 

In Fig. ^] we present the series of graphs for chlo- 
rofluorobenzene equivalent to the fluorobenzene results 
of Fig. O The ionization potentials that correspond to 
the two highest occupied orbitals are 9.37 eV and 9.80 eV, 
meaning that the energy difference between the HOMO 
and the HOMO-1 is 0.43 eV [13 . 

The molecular symmetry is reduced to contain only 
the reflection in the molecular plane as a symmetry op- 
eration {Cs point group) and, consequently, planar nodal 
surfaces can only lie in this plane. However, the two 
active orbitals of chlorofluorobenzene bear some similar- 
ities with their counterparts in fluorobenzene and thus 
our qualitative predictions will also be similar. The sim- 
ilar nodal surface structures can be seen by comparison 
of Figs.[T](c)-(d) with Figs.[T](e)-(f). The HOMO con- 
tains a nodal surface approximately parallel to the xy 
plane and the HOMO-1 contains a nodal surface approxi- 
mately parallel to the xz plane. The similarities between 
the orbitals can be further exploited by comparing the 
asymptotic coefficients which are given in table □ 

When we compare Figs. El and ^| we do indeed find 
alike behaviour. For the second highest orbital the sup- 
pression at X = 0° is now less pronounced. We can ex- 



plain this phenomenon by the fact that the polarization 
is not parallel to a nodal plane but only approximately 
parallel to a nodal surface. 

The most significant difference between the MO-SFA 
LG and the MO-ADK model is found at x = 90°. Here 
the molecular plane lies parallel to the polarization and 
the angle 6 determines how the molecule is rotated in this 
plane. In the MO-SFA LG model the total yield is min- 
imized when 9 > 90° which corresponds to a geometry 
where the chlorine atom is directed along the polariza- 
tion vector. In the MO-ADK model the minimum yield 
is found at 9 < 90° where the fluorine atom is directed 
along the polarization vector. We expect that this dis- 
crepancy arises from the fact that the tunneling theory 
only accounts for ionization in a narrow cylinder along 
the polarization direction of the field, as it is discussed in 
Sec. N~K\ and thus the MO-SFA LG model predicts the 
more accurate results. 



VI. CONCLUSIONS AND OUTLOOK 

In the present work, we have given a detailed account 
of the MO-SFA in the length gauge formulation, and we 
have presented the MO-ADK tunneling theory for non- 
linear polyatomic molecules. We have discussed the va- 
lidity of the present effective single-electron models, and 
we have applied them to molecular systems where we ex- 
pect the approach to be well justified. Indeed a compari- 
son of our theoretical results with experimental ion-signal 
data for ethylene and benzene showed very good agree- 
ment, in particular in the length gauge formulation of the 
molecular strong-field approximation. 

We have presented a detailed study of the characteris- 
tics of orientational-dependent ionization signals for ethy- 
lene, benzene, fluorobenzene, and chlorofluorobenzene, 
and thus covered systems with D2h, D^h, C2v, and Cg 
point group symmetries. Our calculations of angular dif- 
ferential rates for ethylene showed that the distributions 
are largely determined by the nodal plane structure, and 
the comparison of predictions of the MO-ADK and the 
MO-SFA results for the orientational dependent ion sig- 
nal for ethylene illustrated the shortcoming of the for- 
mer model in accounting for situations where the ioniza- 
tion is peaked in other directions than along the polar- 
ization direction. The comparison of the orientational- 
dependent ion signals for the four different molecules 
showed that it is in general the nodal plane surfaces that 
determine the structure in the ionization yields. The 
point group of the molecules, although related, plays a 
less important role. An interesting feature, however, as- 
sociated with the reduction in the point group symme- 
try is the lift of the degeneracy in the HOMO as we go 
from benzene to the substituted benzene molecules. In 
fluorobenzene and chlorofluorobenzene the energy differ- 
ence between the HOMO and the HOMO-1 is compa- 
rably small, and one needs to consider ionization from 
both orbitals. The conclusions for the HOMO-1 with 
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respect to the importance of nodal surfaces remain un- 
changed. If one, on the other hand, speculates about pos- 
sible experimental studies aiming at the separation and 
unique identification of HOMO and HOMO-1 dynamics, 
respectively, we propose to turn to a study of the pho- 
toelectron spectrum or above-threshold-ionization (ATI) 
spectrum. In the fixed-nuclei approximation applied 
throughout this work, we would see an ATI spectrum 
consisting of two series corresponding to electron ener- 
gies which fulfil = nu — Ej,{j) — Up with Ei,{j), 
[j = {HOMO, HOMO-1}) denoting the binding energies 
of the orbitals. Our present results show that depending 
on the molecular orientation, either one of the series will 



lead to the stronger peaks. This effect is of course absent 
in benzene. Beyond the fixed-nuclei approximation, the 
ATI peaks will broaden due to coupling with the rovibra- 
tional degrees of freedom, and the different series may be 
impossible to distinguish. Whether this happens will be 
a study for the future. 
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